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We consider the influence of parton saturation in tlie Color Glass Condensate on 
fSI ' the back-to-back azimuthal correlations of high pT hadrons in pA (or dA) collisions. 

m : 

■ When both near-side and away-side hadrons are detected at mid-rapidity at RHIC 

\ energy, the effects of parton saturation are constrained to transverse momenta below 

Q^' the saturation scale pT < Qs', in this case the back-to-back correlations do not 

Oh: 



X 



disappear but exhibit broadening. However when near-side and away-side hadrons 
are separated by several units of rapidity, quantum evolution effects lead to the 
depletion of back-to-back correlations as a function of rapidity interval between the 



' detected hadrons (at fixed pt)- This applies to both pp and pA (or dA) collisions; 

however, due to the initial conditions provided by the Color Glass Condensate, the 
depletion of the back-to-back correlations is significantly stronger in the pA case. An 
experimental study of this effect would thus help to clarify the origin of the high px 
hadron suppression at forward rapidities observed recently at RHIC. 

Recently, a strong suppression of the high hadron yields has been observed at forward 



as a signature 



rapidities at RHIC lliy,y,y^. Since this effect has been predicted 

of ,„a,.„,„ evolufo. ,„ the Colo. Glas. Condensate flfj y the obse^t.ons have 

excited considerable interest. In this paper we consider an observable which allows to test 
further the origin of the observed effect - the azimuthal back-to-back correlations of high 



Pt hadrons. To do this, we extend the KLMN approach P, [l^ ll^ Q , that has been 
developed to describe the experimental data on hadron multiplicities and the inclusive high 
Pt yields, to the azimuthal correlations 1^ . 

The azimuthal correlations provide a powerful method for the diagnostics of a partonic 
system (for a recent treatment of azimuthal correlations in nuclear collisions, see As 
we will show in this paper, the measurements of the strength of back-to-back correlations 
allow one to tell whether the partonic system under study has reached the density needed 
for the formation of Color Glass Condensate (CGC), or whether it is still in the perturbative 
QCD (pQCD) phase. Indeed, in leading order pQCD a typical hard scattering process at 
high energy is composed of a gluon jet with a large transverse momentum balanced in 
the opposite direction by another gluon jet with transverse momentum {p2,t) which is also 
large and almost compensates the value of pi^t, namely, p2,t — Pi,t = ^1 + ^2*^ \Pi,t\ ( 
see Fig. Q). However, in the CGC phase of QCD the phenomenon of saturation implies a 
different structure of the event: a jet with large transverse momentum can be compensated 
by the production of several gluons with the average transverse momenta which are about 
equal to the saturation scale (Qs) {p2,t ~ Qs)- 
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FIG. 1: Back-to-back correlations 
of produced jets in the perturbative 
phase based on the QCD factoriza- 
tion. 
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In the transitional region dominated by quantum evolution ("Color Quantum Fluid" 
phase or region of extended scaling), the dynamics is driven by the interplay of these two 
mechanisms, which we are now going to discuss in more detail. 

The first one is the production of two gluon jets from one parton shower (see the first 
diagram of Fig. ^ while the second mechanism is the production of two jets from different 
parton showers (see the second diagram of Fig. |21). 

Due to AGK cutting rules [2^, the contribution of the one parton shower production to 



FIG. 2: Two mechanisms of double inclusive production. 



the double inclusive is described by one Mueller diagram of Fig. OJ-a. 
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FIG. 3: Double inclusive cross section for two correlated jet production in the perturbative phase 
of QCD (Fig. 131-a) and single inclusive cross section for a single jet production (Fig. IHl-b). 



The vertex of the gluon emission T^ in Fig. Ol-b (so-called Lipatov vertex) is equal to (see 
for example Ref. Q|) 

2^ 



9 [Pl,t Pl,fJ. 

Pit 



which leads to 



kt {Pi,t - h 
pit 



(1) 



(2) 
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Taking these equations into account one can see that the double inclusive cross section given 
by the diagram of Fig. El-a is equal to (see Fig. El-a) (see Ref . |3| for details) 

1 rfV 



a dyi dy2 d^pi^t (Pp2,t 

2 ^ 2 / ^* 'fprojectileikt, Y - ^i) (ptargeti\Pl,t - P2,t - 4|, 2/2) 



w^i) pipi/ ^\vu-v^A^y-y^^y-) (3) 

It should be stressed that the single inclusive cross section of Fig. Ol-b can be rewritten 
through the same function F^^'^^^ namely 

a dy ^(Ppi " ( aT^"- 1 ) / ^* 'fprojectileikt, Y -y) iptargeti\Pl,t " 2/) 

^4iV.a.^ 1 p,.c.^^^^^^Y - y.y) (4) 



The production of jets from two different parton showers which is described by the second 
diagram in Fig. |21 can be calculated using the Mueller diagram of Fig. |3J It is easy to 
understand that this diagram gives the double inclusive cross section in the factorized form 

1 d'^a 1 da 1 da 

— (01 

a dyi dy2 d'^pi^t d?P2,t cr dyi d^pi^t o dy2 d'^p2,t 
corresponding to the independent (uncorrelated ) production of two gluons with kinematic 
variables {yi,Pi,t) and (2/2,P2,t)- 

Let us now define the correlation function in the azimuthal angle between the two 
gluons as a probability to find a second gluon with rapidity y2 and transverse momentum 
P2^t moving at the angle with respect to the trigger gluon with rapidity yi and transverse 
momentum pi^f Defined this way correlation function has the following form: 

R{(t)]PuP2,yi = 2/2) = (6) 

F'''^H\Pi,t-P2,tlY-yi,y2) + F^^^^(pi,,, F - 7/1, 7/1) F^^^^(p2,., F - 7/2, 2/2) 
/ d(P{F^^CL^\pi,t-p2,t\,Y-yi,y2) + F^^^^(pi,t, 1^ - 2/1, 2/1) F^^^^(p2,t, F - 1/2, 2/2)) ' 
The azimuthal angle dependence originates only from the production of two gluon jets 

from the same parton shower (see Eq. (jS))) while the second mechanism (see the secoond 

diagram in Fig. |2)) leads to the constant background (see Eq. (0)). 
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FIG. 4: The Mueller diagram for the production of two gluon jets from different parton showers. 



The shape of the differential cross sections Eq. Q and Eq. (0) depends crucially on the 
unintegrated gluon densities ip. Here we will use a simple model for these functions adopted 
earlier in Ref. Ili 



5 (1 - xf pt < Qs{x) ; 



(7) 



S 



(1 - xf Pt > Qsix) ; 



In Eq. ((Tj) we neglect, for the time being, the anomalous dimension of the gluon densities 
and use a very simplified assumption about the behavior of reflecting the fact that inside 



the saturation region the density is large and changes slowly j2l|. The numerical factor k 
can be found from RHIC data on ion-ion collisions, but the value of R (see Eq. (jHl)) does not 
depend on it. We introduce, as before |^, 14, 1^, ll^, 12|, the factor (1 — xf which describes 
that the gluon density is power suppressed at x ^ 1 according to the quark counting rules. 
However if we restrict ourselves to calculation of the correlation function at ?/i = ?/2 = 
then the influence of these factors is very small. In general, one can expect this ansatz for 
if of Eq. ((7j) to be a rather crude model but it turns out to be quite successful in describing 
the data on rapidity and transverse momentum distributions j22| . Therefore, we hope that 
our calculations will provide a reasonable guideline for the experimental measurements. 

For the numerical estimate we take the value of saturation scale Ql{A; y = 0) = 2.44 GeV'^ 
for gold at W=200 GeV and Ql{p',y = 0) = 0.26 GeV^ for the proton in accord with our 
multiphcity calculations for deuteron - gold coUisions ^/T,]. The result is plotted in Fig. 
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FIG. 5: Fig. |31-a : Two gluon jet production in one parton shower normalized in a such way that 
the integral over azimuthal angle is equal to 1. The upper curve corresponds to proton-proton 
interaction, the middle one describes the deuteron-gold interaction while the third curve is related 
to gold-gold interaction. The calculation is performed for two jets with pi = A GeV and P2 = 2 GeV 
which corresponds to STAR experiment measurements [2^. Fig. [51-b: Same as in Fig. [SJ-a, but for 
the sum of one and two parton showers. 

One can see that the azimuthal angle distribution has a maximum at = vr which 
corresponds to the jet produced in the opposite direction to the trigger jet. The width 
of the angular distribution is different for pp and dA processes. Introducing the Gaussian 
distribution 

1 (<l,-nf 

R = , e"^^ (8) 



2 7ra 



with the width a, we can see from Fig. El that for proton-proton collision cr = 0.7 while 
for deuteron-gold interaction a = 0.8. Therefore, the difference between the pp and dA 
processes at mid-rapidity appears quite small. However, for gold-gold interactions the 
distribution in Fig. El appears quite different from the Gaussian given by Eq. (jHl). If fitted to 
the Gaussian form, the value of a appears approximately cr = 1.0 -r- 1.2; however this does 
not characterize well the distribution. This result is easy to understand qualitatively, since 
for large pi and p2 (so large that \pi — P2I > Qs{^) ) the integral over k in Eq. Q stems 
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from the region where k ^ Qs{p) < \Pi ^ P2\ and 



F^''^\\p\,-p,,lY-y„y,) oc ^^4^ - 

{Pi -P2) + 2pi P2 cos<p 1 + P2 

In this formula the typical width of the azimuthal angle distribution is equal to ctypicai = 

{pi — P2Y /'^PiP2- For the kinematics used in the STAR experiment [2^, pi = AGeV and 

P2 = 2GeV, we obtain atypical ~ 0.5 from this simple formula. As was mentioned above, 

the fact that Fig. El leads to a large value of atypical is due to the region of integration over 

In Fig. ini-b we plot R{4>; Pi,P2, Vi = I/2) for = 4 GeV and P2 = 2 GeV as a function of 0. 
R has a meaning of the probability to find a particle with momentum p2 at a definite angle 
if a trigger particle has a momentum pi. The curves in this figure are quite different from 
Fig. El due to the background from the production of gluon jets from two parton showers. 
This background is proportional to N^^^.^ where Npart is the number of participants (see Refs. 

Q| ). In fact, this background becomes so large for gold-gold collision that we cannot 
see the azimuthal correlations. However, for proton - proton and deutron-gold collisions the 
background is not so large and we see a clear peak at = vr. 

However, experimentalists in their analysis sometimes substract the fiat background - the 



"pedestal" 2J]. Doing such a subtraction for Fig.|31-b and normalizing the - distribution to 
the unit area we arrive at the correlation functions shown in Fig. IHl It should be stressed 
that the ratio of the areas is equal to pp/pA= 1.16. Thus one can see that we get quite 
similar - distributions for proton-proton collisions and for deuteron-gold ones, in (at least 
qualitative) agreement with the data jzJ ]. 

To make the main features of our approach even more transparent, let us use the fol- 
lowing simple model for the unintegrated structure function ip which allows us to make all 
calculations analytically: 



2 



One can see that for pt <^ Qs and for pt ^ Qs we have the function cp of Eq. ^ , while in 
the region of pt ^ Qs f of Eq. (fTT1|) is different. The advantage of Eq. (|Tn|l is the fact that 
we can calculate F^^'"^ of Eq. Q and Eq. Q analytically. Indeed, introducing Feynman 
parameters and taking integral over kt in Eq. Q and Eq. (0)) we obtain 

i^'^''^(Pt,Q.,g.) = vr [' dt , ^ -— ^ (11) 

Jo t{l-t)pi + t{Ql-ql) + g2 
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FIG. 6: The azimuthal angle correlations after substraction of the background. The upper curve 
corresponds to proton-proton interaction while the lower one describes the deuteron-gold interac- 
tion. The calculation is performed for two jets with pi = A GeV and p2 = 2 GeV which corresponds 

Q 



to STAR experiment measurements 



TT ^^f Pt + Q'i + ql + h{pt, Qs, Qs 



2 h{pt, Qs, qs \Pt + Ql + ql - h{pt, Qs, 
where Qs and qs are saturation momenta of projectile and target and 



h{pt,Qs,qs) = ^{pl + {Qs-qs?{pl + {Qs + qs? (12) 
One can see that the two scales governing the azimuthal angle dependence emerge: 



(-) _ \Fi — P2) -r — qs) (+) _ KPi — y2) t yy^ -r gs) 

In the case of ion-ion collisions the value of cr typical is rather large, and this leads to a wide 
distribution in 0. Using this simple model, we recalculate the results presented in Fig. 
Fig. |S]-b and Fig. IHl using Eq. (fTTJ); the results are presented in Fig. Q Fig. [7|-b and Fig. |H1 
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Fig. Qb 



FIG. 7: Fig. [3-a: Two gluon jet production calculated in the simple model (see Eq. for the 

unintegrated stricture function (/? given by Eq. ()1U() in one parton shower normalized in a such way 
that the integral over azimuthal angle is equal to 1. The upper curve corresponds to proton-proton 
interaction, the middle one describes the deuteron-gold interaction while the third curve is related 
to gold-gold interaction. The calculation is performed for two jets with pi = A GeV and p2 = 2 GeV 
which corresponds to STAR experiment measurements 2J]. Fig. [TJ-b: Two gluon jet production, 



calculated using Eq. in one and in two parton showers normalized in a such way that the 

integral over azimuthal angle is equal to 1. 

Examination of these results leads us to the conclusion that our approach explains, at least 
on a qualitative level, two of the most striking experimental facts: (i) the strong suppression 
of the azimuthal back-to-back correlations in gold-gold collisions and (ii) the close similarity 
of these correlations for the proton-proton and deuteron-gold interactions. (The final-state 
interactions of the jets in ion-ion collisions are expected to broaden the observed azimuthal 
correlations even more.) 

It is interesting to see if we can reproduce the measured widths of the 0-diastribution 
quantitatively. As one can see from Fig. IHl the calculated width is larger than the experi- 
mental one. One obvious reason which can contribute to this difference is the fragmentation 
of the jets: indeed, we have computed the correlation function for gluons, not for the mea- 
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FIG. 8: Fig. IHl-a: The azimuthal angle correlations, calculated using Eq. (|1H) . after subtraction 
of the background. The upper curve corresponds to proton-proton interaction while the lower 
one describes the deuteron-gold interaction. The calculation is performed for two jets with pi = 
4 GeV and p2 = 2 GeV which corresponds to STAR experiment measurements Q Fig. Hb: The 
azimuthal angle correlations of produced hadrons after substraction of the background. The upper 
curve corresponds to proton-proton interaction while the lower one describes the deuteron-gold 



interaction. The calculation is performed for two 
corresponds to STAR experiment measurements 2 
of Eq. Q with a = 04. Fragmentation is included as described in the text 



ets with pi = 4 GeV and p2 = 2 GeV which 
. The dotted line is the Gaussian distribution 



sured hadrons. As a result of fragmentation, the fall-off of the inclusive cross section (and 
also of p^^(^^ in Eq. Q and Eq. (0)) with the transverse momentum becomes more steep. 
To take account of this, we use the fragmentation function of gluon jet to hadrons from Ref. 
[2^ ■ One can see that the resulting width of the distribution in Fig. |Hlb is indeed much 
smaller than in Fig. IHl and is now close to the experimental one. One can also see that the 
real distribution is not quite Gaussian. 

One of the possibilities to check the validity of our approach is to measure the az- 
imuthal correlations in the kinematics when both of the detected hadrons are measured 
at forward rapidity. Since the nuclear saturation momentum increases with rapidity. 
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FIG. 9: The azimuthal angle correlations of produced jets at forward direction with rapidity 
rj = 3.8. The value of the trigger momentum is taken 1.5 GeV while the hadrons produced in the 
backward direction were integrated over transverse momenta from 0.2 GeV to 1.5 GeV. 



Qs{A,y) = Qs{A,y = 0) exp[A?/] with A ~ 0.25 -r 0.3 [2^, we expect a wider distri- 
bution in the azimuthal angle, than at y = 0. Fig. El shows our prediction for deuteron-gold 
and proton - proton scattering for y = 3.8. 

A most interesting opportunity to investigate the CGC dynamics in the quantum domain 
is to study reactions 27[ in which the trigger hadron with transverse momentum pi is 
generated in the forward direction ( say, yi = 3.8) while the recoiling particle(s) produced 
at the azimuthal angle is detected in the central rapidity region (?/2 = 0, see Fig. irUll . 
The interest in this kinematics stems from the fact that a large rapidity interval Ay = 
Vi — 2/2 between the detected particles enhances the effects of quantum evolution, since the 
probability of gluon emission is proportional to Ay. In fact, the study of gluon jets 
separated in rapidity has been proposed by Mueller and Navelet |2S,] as a way to investigate 
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the properties of BFKL evolution. In the case of nuclear target, the quantum evolution 
enhances the influence of the saturation and extends it to larger transverse momenta. 
The double inclusive cross section for such kinematics can be written in the form: 

1 da ^ f 4.N,as \^ 1 1 

a dyidy2 d^pi^t d'^Pt,2 -iJ pit p\t 

d^ kt , d^ k[lf projectile (^t , Y " Vl) f^^^^ {Vl - 1/2, h, K; 0) !f target {K, ^2) 

where the BFKL scattering amplitude for two gluons. This amplitude takes 

into account the emission of gluons with rapidities between yi and 2/2 (see Fig. irUj) . This 
amplitude can be written as series 

00 

^^^^^ {yi - y2, kt, K; 0) = E Cos{n 0) ^™ (^2 - yu h, K) (15) 

n=0 

where (f) is the azimuthal angle and ip^^^^ are the eigenf unctions of the BFKL equation, 
which at high energies behave as 

^™(2/2-yi,A;*,A;0 ^ e^" (^-^'^) (16) 

where 

^„ = ^(V'(i) - + (17) 

However, in Eq. (jl5p only the first term has a positive intercept {ujq > 0) while all other 
terms fall off as a function of energy. Indeed, Eq. (fTTj) gives 

4 \n2 asNc _ 1.22 asN^ 2.56 asN^ 

Uq — ; Ui — ; U2 



TT TT vr 

We thus replace (p^^^^^ by sum of two first terms with n = and n = 1, since other terms 
are supressed at large value of \yi — y2\. 

Fig, im shows the normalized azimuthal distribution of the jets when the recoiling jet has 
?72 = and 0.2 GeV < P2 < 1.5 GeV while the trigger jet has rji = 3.8 and pi = 1.5 GeV 
assuming that both of the jets are produced from one parton shower (see Fig.|31-a and Fig.^J. 
One can see that in both cases (for proton - proton and deuteron-gold collisions ) we expect 
sufficiently strong correlations. For deuteron-gold collision the width of the distribution in 
the azimuthal angle is only 30% larger than for the proton -proton scattering. 

This result is not final though, since the main difference between the two cases is in the 
independent production from two parton showers ( see Fig. HI). In the CGC phase we have 
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FIG. 10: The BFKL emission in the 
double inclusive cross section. 



a saturation of the parton densities which is expressed in our assumption for the functions 
<y9's (see Eq. (0)). In this region the diagrams for independent two jet production (see Fig. |3] 
) are much more important than the production of two jets from the single parton shower. 
Indeed, diagrams of Fig. OJ-a and Fig. lead to the cross section of the order of 1 while the 
independent production given by Fig.Hlleads to the cross section of the order of {1/ asiQs))'^ ■ 
In reality, Fig. ITTlb shows that the azimuthal correlation in deuteron - gold collision leads to 
a maximum around (p = tt which the height of only 17%. For proton -proton we still expect 
a sizable effect of around 57%. One can see therefore that the quantum evolution effects in 
the CGC indeed induce a large difference between the back-to-back correlations expected in 
pp and dA collisions. 

Let us now discuss the main uncertainties involved in our calculations. Our estimates 
of the independent production in comparison with the production in one parton shower 
have large errors because they involve the normalization of the inclusive cross section. The 
relative contribution of the diagram of Fig. 0] to the diagram of Fig. Ol-a (we define Z as 
a numerical factor which should stand in front of the diagram of Fig. 0] which is defined 
as the square of expression given by Eq. Q.) can be expressed through the normalization 
constant that we have introduced earlier in describing the multiplicity in d A collisions 17 1. 
Denoting this constant by C we have for the normalization of the relative contribution (Z) 
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FIG. 11: Fig. Illl -a: The azimuthal angle correlations of produced jets in one parton shower at 
forward direction with rapidity t/2 = and with 1.5 GeV > P2 > 0.4 GeV when the trigger jet is at 
iji = 3.8 and with transverse momentum pi = 1.5 GeV. Fig. llll -b: The azimuthal angle correlations 
of produced jets at forward direction with rapidity r/2 = and with 1.5 GeV > P2 > 0.4:GeV 
when the trigger jet is at 771 =3.8 and with transverse momentum pi = 1.5 GeV . Both, production 
from one parton shower and from two parton showers are taken into account. 



the following expression 



Z 



(18) 



9asiQs) rijet 

where rijet is the hadron multiplicity of jet with transverse momentum Qg. The uncertainty 
in calculation of Ujet is large and in our numerical estimates we take rijet = 1-5. For proton 
-proton collisions we introduce an additional normalization factor. We need it to describe 
the relation beetwen the parton density and the saturation momentum for this case since we 
cannot trust the geometrical estimates of the area of interaction in this discussed in 

Q. 

It is important to note that the uncorrelated production can be subtracted from the 
data experimentally since the inclusive cross section has been measured. Therefore, the 
correlation can be attributed to the diagrams of Fig. ^1 In Fig. IT^a and Fig. IT^b we 
illustrate the azimuthal angle distribution for two different kinematic regions. 
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Fig. [T2Ja Fig. [T^Jb 

FIG. 12: Fig. I12t -a: The azimuthal angle correlations of produced jets at forward direction with 
rapidity r]2 = and with p2 = 1.5 GeV when the trigger jet is at r/i = 3.8 and with transverse 
momentum pi = 1.5 GeV. Only production from one parton shower is taken into account (see 
Fig. llUj) . Fig. 1121 -b: The azimuthal angle correlations of produced jets at forward direction with 
rapidity r/2 = and with p2 = 2 GeV when the trigger jet is at r/i = 3.8 and with transverse 
momentum pi = 2 GeV. Only production from one parton shower is taken into account (see 
Fig.Hni). 

These figures (Fig. [T21a and Fig. fT^b ) should be compared with Fig. [T21 which gives 
the azimuthal correlations when both triggers have the same rapidity {rji = r]2 = 3.8) and 
originate from the same parton shower. 

In fact one can substract the independent production of two hadrons just by measuring 
the inclusive cross section. The remaining emission in one parton shower then depends both 
on the saturation momenta for colliding hadron and/or nuclei and on the BFKL emission 
in the rapidity interval rji — 772- The prediction for this process is theoretically reliable and 
could provide the information on the values of the saturation momenta. 

To summarize, we have found that hadron azimuthal correlations provide a stringent tests 
of the parton saturation in the Color Glass Condensate. When both hadrons are detected 
at central rapidity, our results show that pp and dA correlations are quite similar, with some 
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FIG. 13: The azimuthal angle correlations of produced jets at forward direction with rapidity 
r/2 = and with pi = p2 = 2 GeV when the trigger jet is at rji = rj2 = 3.8. Only production 
from one parton shower is taken into account (see Fig. llUj) . 

broadening in the latter case, whereas the correlations in AA are suppressed, mainly due to 
the large independent production background. A most interesting possibility is provided by 
the kinematics in which one of the hadrons is detected at forward rapidity and another at 
central rapidity. In this case the effects of quantum evolution in the CGC on the correlation 
function are strong, and enhance the influence of the saturation boundary. This leads to a 
significant difference between pp and dA correlation functions. 
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